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Abstract 

In this work we study numerically a short range p-spin glass model in three dimensions. 
The behaviour of the model appears to be remarkably different from mean field pre- 
dictions. In fact it shares some features typical of models with full replica-symmetry 
breaking (FRSB). Nevertheless, we believe that the transition that we study is intrin- 
sically different from the FRSB and basically due to non-perturbative contributions. 
We study both the statics and the dynamics of the system which seem to confirm our 
conjectures. 



PACS Numbers 05-7510N 



1 Introduction 

One of the presently unsolved theoretical issues concerning spin glasses is to obtain a satisfying 
description of these systems in finite dimension. By now, a quite complete picture has been 
achieved for spin glass models at mean field level [0. Quite generally, one finds a critical 
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temperature T c which separates an ergodic phase from a low temperature spin glass phase in 
which ergodicity is broken. Below T c the spins are frozen in their local orientation, giving rise, 
because of the disordered nature of the interactions, to no total spontaneous magnetization. 

This transition is well described by a spontaneous replica symmetry breaking scheme 
In the replica framework the order parameter of the spin glass transition is a function q(x) 
defined on the real interval [0, 1] that changes from a high temperature constant form q(x) = go 
to a low temperature function that can assume different values. The physical meaning of the 
function q(x) is the following: at low enough temperatures the free energy landscape has many 
local minima separated by free energy barriers that grow, asymptotically, as the size of the 
system. Below the critical temperature the system breaks the ergodicity and remains confined 
into one of these valleys. The values assumed by the order parameter function q(x) are the 
possible values of the overlap q aT = N~ l Y.f &i T i among two copies (replicas) of the system 
which are independently equilibrated in the same free energy landscape. 

If the two replicas freeze into the same pure state then the overlap will be the maximum 
allowed {qEA) and it will provide a measure of the width of the state. The size of the interval 
where q(x) = q gives the probability that the two replicas are found with an overlap q. More 
precisely, one can define a quantity P{q) = where x(q) is the inverse function of q(x). The 
function P(q) is the probability distribution of the values of the overlap and it contains the 
same amount of information as the q(x). At high temperatures one has a constant q(x) = qo 
which means that the P{q) is a delta function on go- 
Two main classes of spin glass models can be distinguished depending on the character of 
their phase transition. 

• A first class of models undergoes, in mean field, to a full replica symmetry breaking 
(FRSB) transition in which the function q(x) changes continuously with the temperature. 
The low temperature solution is a continuous function in which all values of the overlap 
between g and qEA are attained. Correspondingly, the P(q) will change to a non-zero 
continuous function between g and qEA with an additional delta function on q E A- 

This transition is continuous in the sense that qEA — Qo is small for T close to T 6 7. 

These models describe with success the physics of real spin glasses. This kind of transi- 
tion is second order. At a critical temperature T c the spin glass susceptibility diverges 
and long range correlations establish. 

• A second class of spin glasses has a discontinuous low temperature function q(x) that 
can assume a finite number of values. Usually q(x) assumes only two values go an d 
gi(= qEA) an d the transition involves one step of replica symmetry breaking (1RSB). In 
most cases the transition is discontinuous in the sense that g — qEA — 0(1) as soon as 
T<T C . 

The low temperature P(q) can be written in the following form 

P(q) = mS(q ) + (1 - m)5(q 1 ) 

where m is a parameter which is smaller than unity and it decreases with the tempera- 
ture. 
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In this second class of models the distribution P(q) is zero between qo and q\ which are 
the only allowed values for the overlap. This is very different from the case of models 
of the first class where the low-temperature P(q) shows a plateau from zero to qEA and 
there is no 5(qo). In absence of magnetic field one always has go — 0. We shall stick to 
this case in the following. 



Another important feature of mean field spin glasses is that the dynamical equations for 
the two-time correlation functions show a singularity at a critical temperature Td > T c . It 
turns out that for models of the first class the dynamical and statical critical temperatures 
coincide, while in models of the second class one has > T c . The existence of a dynamical 
singularity at a higher temperature than the static critical temperature is an aspect of these 
models which is shared by structural glasses 0. 

The mode-coupling (MC) theory describes quite well the dynamics of supercooled liquids 
[0 approaching their structural glass transition. The MC dynamical equations predict a 
singularity at a temperature T mc which is higher than the glass temperature T g predicted for 
real glasses in a static Adam-Gibbs scenario. 

There has been recent effort to approach the structural glass transition by means of the 
techniques developed for disordered systems such as spin glasses j|, |[]. An effective disorder 
is in fact induced in structural glasses by the slow dynamics of the microscopic variables. 

So, despite the absence of quenched disorder in the Hamiltonian and the possible presence 
of a crystalline groundstate (which however corresponds to a golf-shaped minimum in the free- 
energy landscape and can therefore been considered irrelevant to the physics of the system), 
it is believed that the behaviour of structural glasses is very similar to the one of 1RSB mean 
field spin glass models. Recent numerical results seem to support this hypothesis [||- [pTOf] . Of 
course, in finite dimensional glasses, the effect of the decay of metastable states by activated 
processes have to be taken into account, and are an element of difference with respect to mean 
field spin glasses. An interesting issue is to generalize long range models of spin glasses with 
1RSB to models with short range interactions to see if the finite dimensional effects are similar 
to the case of real glasses or not. 

In this work we present a numerical study of a 3d short-ranged version of the p-spin model 



JL1], [12|, [13|, [14|. The model is provided of a parameter M that ensures the mean field solution to 
be exact when M is sufficiently large. A perturbative analysis around the mean field solution 
predicts a behaviour qualitatively similar to the MF ||11|| . Nevertheless, numerical simulations 
give evidence for a diverging correlation length when approaching the critical temperature 
from above [|l^ . 

We study the model in d — 3 for p = 4 and M — 3. We look at the behaviour of the 
function P(q) which appears definitely non trivial at low temperatures. The form of the P{q) 
that we find seems to be relevantly different from the shape that one has in the case of the mean 
field p-spin model. Our results indicate the occurring of a continuous transition. Furthermore, 
a numerical study of the dynamics of this model seems to provide evidence for a dynamical 
transition at the same critical temperature of the static one as in the cases of models with 
a continuous transition. Nevertheless we believe that this transition is very peculiar and is 
not due to an instability of the saddle point as in the case of models with FRSB. We address 



3 



the peculiar nature of this transition to the formations of large regions of highly correlated 
spins which are localized by the quenched disorder and whose typical size grows as T — > T c . 
Before presenting the short range model and our numerical results, we will devote the following 
section to describing the known results for the 1RSB class of mean field spin glasses and to 
do so we will focus on a particularly simple model belonging to this class. 



2 The long-range £>-spin glass model 

In this section we shall briefly resume the known results for a long range model that has been 
extensively studied in the past and, because of its simplicity, is a recommendable paradigm 
for mean field spin glasses undergoing a 1RSB transition. 

The model in zero magnetic field is defined by the p-spin interaction Hamiltonian 



H p ({a}) = - 



(l<i 1 <i 2 <---<i P <N) 



(1) 



The spin variables can be either Ising spins (cr, = ±1) or continuous soft spin variables 
which have to be constrained by some normalizing condition such as 



N 



i=l 



N, 



(2) 



where the interactions Ji 1: i 2 ,---,i p are random variables distributed with 





2 


irJ 2 p\ 


exp 



(■Til,' 



1 N P- 



J 2 P \ 



(3) 



The scaling of the variance with N p ~ l ensures the free energy to be extensive. The mean 
field solution of the model is exact when the limit N — * oo is taken since the interactions are 
long ranged. The physics of the Ising and of the spherical p-spin glass model is essentially the 
same if p > 2. Here we will review some well established results in the spherical case which 
can be exactly solved. For p = 2 this model does not belong to the 1RSB class, and exhibits 
a continuous transition. 

The model presents a 1RSB discontinuous transition at a temperature T c . One can dis- 
tinguish two other relevant temperatures T d and T m , where 

T c <T d < T m . 

The physics of the model is characterized by the appearance of a large number of solutions 
different from the paramagnetic one. We can schematically resume the behaviour of the model 
with respect to temperature as follows 



For T > T rn the model is in a paramagnetic phase which is the only minimum of the free 
energy. 
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• For Td < T < T m the free energy landscape changes. The physics of the model is still 
dominated by the paramagnetic state, but other local minima of the free energy appear. 
These minima have a lower energy than the paramagnetic state but their free energy is 
higher and therefore they are still irrelevant in the thermodynamic limit. There is an 
exponential number of these minima. One can define the complexity £ as the logarithm 
of the number of states that becomes extensive below T m . 

The free energy of the ensemble of these states is therefore 

F S = E S -TS- TS (4) 



Above Td one has F s > F p where F p is the free energy of the paramagnetic state. So in 
this regime the dominant solution is still the paramagnetic one. 

• In the whole range of temperatures such that T c < T < T d one finds that the free energy 
of the paramagnetic state equals the free energy of the ensemble of local minima i.e. 

T = F 

It turns out that at the temperature Td the solutions to the dynamical mean field equa- 
tions change. In this sense Td is the dynamical critical temperature. This corresponds 
to the fact that, if the system starts from a random initial configuration, the system will 
eventually fall into a metastable state where it will be trapped for an infinite time. 

• For T close to T c , E(T) vanishes linearly with T — T c . At T = T c the number of local 
minima becomes non-extensive and a static transition of the 1RSB kind occurs. 

As the system is cooled below T d a dynamical transition occurs. One can easily write some 
closed equation for two time quantities such as 



t w , t w 

If the system has reached equilibrium, the two time quantities depend only on the differ- 
ences between the two times and are related by the equation 

GO) - (7) 

These two properties, time-translation invariance (TTI) and the fluctuation-dissipation 
theorem (FDT) are properties of equilibrium dynamics. 

Below Td the system has a regime of slow dynamics where, for large time differences t ~ t w , 
the relaxation of C and G is slower the larger t w is i. e. the older the system is. In this regime 
one can write the two-time quantities under the form 



1 

i 

1 " d(<Ti(t + t w )) 
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C(4 + t „, U=(5 p!±«) G(t + t ^ ) = T l_G(^±f), (8) 

Relation (^) is not valid under T^, but a generalized FDT can be written introducing a 
function X tw (C) ||5j by the relation 

G(() = -^^). (9) 

For large waiting times one has that X tw (C) — > X(C). The function X(C) generally 
characterizes the type of aging dynamics of the model. In models where Td = T c (i.e p = 2) 
one has that X(C) = x(q) while this equation does not hold in general though the two 
functions remain equal in structure ||16|| . The dynamical transition in mean field spin glasses 
is interpreted as due to the freezing of the system into one among the O(expiV) metastable 
states, which for a long range model have an infinite lifetime in the iV — > oo limit. 

The presence of metastable states with infinite time life is just an artifact of the mean field 
approximation and we do not expect Td to mark an effective dynamical transition in a finite 



dimensional model. Correspondingly one should find [18| that the dynamical function X tw (C) 
becomes equal to the static x(q) = ^dq'P^q') in the limit of t, t w — > oo (in which C — > q) 
also for 1RSB models. Nevertheless, it should be still possible to observe a Td higher then the 
critical temperature T c . In a finite-dimensional system that has a 1RSB mean field solution 
Td should mark the onset of a two step relaxation process: the mean field-like relaxation of 
the system within a metastable state and the following decay of the metastable state due to 
activated processes. 



3 The model 

The model that we study in this work is a short range version of the Ising p-spin glass. 
We believe that for our purposes the Ising and spherical version of the p-spin are essentially 
equivalent. 

The model is defined as follows. We consider a cubic tridimensional lattice of side L. 
On each site of the lattice we put M Ising spins. The total number of spins is therefore 
N = M L d . Each spin interacts with M — 1 spins on its same site and with 2dM spins on its 
nearest-neighbour sites. The interactions J^Zil are quenched random variables. 

The Hamiltonian of the model can be written as 

L d M 

n P ({*})= E E JtX^ 1 ( 10 ) 

<ii,---,i p > Zi,---,Zp=l 

By 53<ti,...,i p > we m ean the sum over all the sites of the lattice taking, for each couple of 
adjacent sites i and j, p — k of the ii, ■ ■ ■ , i p indexes equal to i and k indexes equal to j. In 
other words, for each nearest neighbour sites i and j, every interaction involves p — k spins of 
site i and k spins of site j with k running from zero to p. 
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We consider discrete (±1) spin variables and we call <Ji r lr the I 1 / 1 spin of site i r with l r 
running from 1 to M. 

For large M, each spin interacts with a large number of nearest neighbours and, in the 
limit M — > oo, the mean field approximation has to be exact. An analytical study of the 
model for large M has been performed in [I1J where a Gaussian P(Jil"'/) is considered. The 



mean field solution of the model coincides with the solution of the long range Ising p-spin. 

One could hope to peep into the finite-M case by perturbing around the M — > oo limit. 
Calculating the 0(1/ M) corrections to the free energy and the Gaussian propagators one sees 
that a perturbative approach shows a transition which is very similar to the MF one, with a 
discontinuous q(x) and no diverging low-momentum static propagators. In the model that we 
simulated we chose the following distribution for the couplings 

P(J) = ^,1 + ^-1- (11) 

The fact that the distribution of the couplings is non-Gaussian should not be source of 
relevant differences from the analytical calculation. Any result in disagreement with the results 



of [1 i] should be imputed to finite dimensional effects which can not be reached with 0(1/ M) 
calculations. 

Here we take M = 3 and p=4. Since p is even, the Hamiltonian is invariant for inversion of 
all the spins ( {<?/} — > {— CTi 1 } ) and the P(q) is correspondingly symmetric, i.e. P(q) = P(—q). 

As usual in spin glass simulations we consider two replicas (with the same configuration 
of disorder) of the system, with spin {a} and {r} respectively, that evolve simultaneously and 
independently. We define 

-i L 3 i L 3 1 M 

q = ±Yq=±y±Y <t i t j (12) 

i=i i=i i=i 

The order parameter probability distribution P(q) is then given by 



P(q) = (S(q-Q)) (13) 

where (•) means thermal average and (■) means average over disorder. 

When looking at the out of equilibrium behaviour of the system it is useful to consider 



17H the staggered magnetization 



1 L 3 -I M 

^s[h}(t) = -Y,(-j:a l t (t)), (14) 



that is related to the response function G(ti,t 2 ) by 



G( ^H« (15) 



One has therefore 



m s [h\(h) = f 1 dt J^M = f 1 dt 2 G(tt, t 2 )h(t 2 ) + 0(h 2 ), (16) 
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which is just the linear-response theorem. 

By applying the generalized FDT relation (§) one obtains: 

m[fc](ti) ~ p f dt'X[C{t u t 2 )] dCitlM) h{t 2 ). (17) 

J~ oo (712 

If the small perturbing magnetic field is turned on at the time t w to a constant time- 
independent value ho we get 

m a [h]{*i) * 0h o J^d1fX[C(t ll t 2 )] dC( & ,t2) = 

(18) 

= Ph Sc(tutv,) duX M = P h oS(C), 

where we have used C(t\, t^)—! for t\ — t 2 . 

In the case of finite dimensional models, for large t w , t\ = t + t w , one has fl8| C — > 
q, X(C) — > x(q), where x(q) = J q dq'P(q'). This means that the staggered magnetization 
measured by out of equilibrium simulations is related to static quantities and more precisely 

S(C)^y(q)= ( l dq'x{q') (19) 



We are therefore allowed to confirm results on the equilibrium probability distribution of the 
overlaps P(q) by looking at the dynamical behaviour of large systems. 



4 On the behaviour at the equilibrium 

4.1 Simulations 

In order to measure equilibrium quantities we use Parallel Tempering (PT) [0. We consider 
size ranging from L=3 up to L=6 (iV=648), simulating contemporaneously two independent 
sets of n replicas that move between f3 min = 0.14 and (3 max = (3 min + nA/3=0.5 (i.e. down to 
T min =2.0), where A/3=0.4 for L=3, 4 and Af3= 0.2 for L=5, 6. We perform 2 18 PT steps for 
L=3, 2 19 for L=4 and 2 21 for L=5, 6, that means up to 2 millions of MC steps for each of the 
2 n=38 replicas. 

Thermalization is checked in different ways: 

• We measure all the relevant quantities during the whole run (in each 2* - 2 t+1 interval), 
looking at the convergence to equilibrium and checking that there are no shifts of the 
mean values in the last intervals. In particular, we do not find evident changes in the 
behaviour of P(q). The presented data are the ones collected in the last half of the run. 

• We check that each replica moves several times from an extrema of the temperature 
range to the other and back. 

• We evaluate the specific heath c both using c = d < e > jdT and using T 2 c =< e 2 > 
— < e > 2 , checking the compatibility of the results (see [Fig. |l|). 

Errors are estimated from the fluctuations between different disorder realizations. We 
consider 400 samples for L=3, 300 for L—4, 240 for L=5 and more than 100 for L=6. 
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4.2 Results on xsg an d P(q) 



We present in [Fig. 1] data on the energy and the specific heath as a function of T for the 
different sizes considered. Finite size effects are quite evident for L=3 while the observed 
behaviours result very similar for larger L values. 




Figure 1: Data on the energy density e (left) and on the specific heath c (right) as a function 
of T for L=3(+), 4(x), 5(*) and 6(D). Lines are only to join neighbouring points. 

The spin glass susceptibility xsg = M L 3 (q 2 ) is plotted in [Fig. 2] for the different sizes 



considered. The obtained behaviour confirm [O] the presence of a non-zero T c thermodynam 



ical transition which appears to be of second order. Data agree well [Fig. 3] with the finite 
size scaling law 



XSG (L, T)/L 2 ^ = xsg {L l, \T - T c )) , (20) 
where we use T c ~ 2.6, r\ ~ and v ~ 1. The value of v results slightly higher than the 



previously estimated [12 v ~ 2/3, this being probably due to some corrections to finite size 
scaling. 

In [Fig. 4] we plot P(q) for L=5 at different temperatures, from T=3.33 down to T=2.0. 
The behaviour emerging from these data results quite puzzling. 

When approaching T c , P(q) becomes clearly non Gaussian but there is no evidence for the 
transition being discontinuous in the order parameter. Our data are compatible both with 
qi=0{T — T c ) at T<T C and with the possibility of a finite q%=0(l), hidden by finite size effect. 

At lower temperatures the situation appears still less clear: From the mean field theory 
we expect only two possible values of the overlap, i.e. go — and q% which increases when 
lowering T, whereas we find a P(q) which is clearly greater than zero in all the range [—51, q\] 
also at the lowest temperature considered T=2.0. Moreover we observe, as expected, the 
weight of the peak in q to decrease when lowering T but it widens, although never disappear 
completely. 

We could suppose that these are finite size effects. It should also be noted that in the mean 



field spherical p-spin all the values of the overlap are allowed within metastable states [|20] and 
the contribute of these states is expected to be negligible only at very large N. Nevertheless 
these effects do not appear to decrease as L increases. 
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Figure 2: Data on the spin glass susceptibility xsg as a function of T for L=3(+), 4(x), 5(*) 
and 6(d). Lines are only to join neighbouring points. 
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Figure 3: Data on the spin glass susceptibility xsg/L 2 as a function of (T — 2.6) L for L=3(+), 
4(x), 5(*) and 6(D). 
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Figure 6: Data on P(q) for L=6 (118 samples) at the lowest temperature considered T = 2.0 

as obtained in the last half of the run (right) and in the previous interval (i.e. in the previous 
2 20 _ 2 i9 PT steps ^ 



As it can be seen from [Fig. 5], where we present data on P(q) at T=2.0 for the considered 
sizes, the minimum value P m in = mrn ge[-gi,(?i] P(o) f° r L=3 (N=81) is compatible within the 
errors with the one for L=6 (iV=648). Moreover the broadening of the peak in g when 
considering larger sizes is very evident if we compare data for L=3 with the ones for L=5. 

The comparison between the L=5 and L=6 data on P(q) results puzzling since the weight 
of the peak in go seems to be larger in the last case. On the other hand it is well known that 
a correct estimation of -P(O) in spin glass simulations requires great care since the possible 
presence of non well thermalized samples causes the overestimation of this value. 

In order to clarify this point we plot in [Fig. 6] P(q) for L=6 at T=2.0 in the last half 
of the run (the same data that in [Fig. 5]) and in the previous interval (i.e. in the previous 
2 20 — 2 19 PT steps). Although data are compatible within the errors, the shift of the peak in 
go is quite evident (this is an effect that we do not observe in L=5 data) and it shows that 
some of the samples are still not perfectly equilibrated at the beginning of the last half of the 
run. We are therefore suspicious that 2 23 PT steps would be needed to be sure of achieving a 
complete thermalization for L=6 down to T=2.0, which means unfortunately too much CPU 
time. 

If the peak in go would disappear completely when considering larger size, the P(q) would 
display a FRSB-like behaviour at low temperatures, i.e. two peaks in ±qi separated by a 
continuous plateau. Nevertheless, by looking carefully both to L=5 and L=6 data at low 
temperatures and to the previous [Fig. 5] , we note that the broadening of this peak is related 
to the appearance of smaller, not well pronounced peaks near go and that the behaviour of 
the -P(g) near and below T c results very different from the one usually encountered in short 
range models with FRSB [2jfl . 

To emphasize the last observation, we present in [Fig. 7] data on the cumulant 




(21) 



This is the well known dimensionless parameter usually successful for evaluating T c , since g(T) 
should go to zero above T c in the thermodynamic limit and the curves for different values of 
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Figure 7: Data on g(T) for L=5 (240 samples) at different temperatures. From left to right 
and top to bottom T=3.33, 2.94, 2.63, 2.38, 2.17 and 2.0. 



./V should cross at the critical temperature because of finite size scaling. In the FRSB mean 
field models as well as in their short range versions (EA models) g(T) for finite N is an 
always positive function of T that increases quite regularly when lowering the temperature. 
The behaviour displayed in [Fig. 7] results therefore in evident disagreement with the one 
expected in the usual FRSB case: we find a g(T) which is a non monotonic function of T, 
negative on a large part of the range of temperatures considered. This accounts for the fact that 
the function P(q) develops non-Gaussian tails in a range of temperatures above T c detecting 
the presence of non perturbative effects (see the related discussion in the conclusions). The 
shown data correspond to L=5 but the behaviour is qualitatively similar for the other L values 
that we have considered and it suggests that this model belongs to a different universality class 
from the EA spin glass one. 

In conclusion, more extensive simulations are necessary to achieve a better understanding 
of the short range p-spin glass. Nevertheless it seems quite evident from our data that we 
are looking at a model characterized by a P(q) that is in disagreement with the one of 1RSB 
mean field models but that appears also very different from the one usually encountered in 
FRSB models. This is a quite new result that may have interesting physical interpretations 
as we will discuss in the following. 

5 On the dynamical behaviour 
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5.1 The approach to the equilibrium 



We have already pointed out that the infinite time life of metastable states is an artifact of 
mean field theory. We therefore do not expect to find a dynamical transition at Td > T c in a 
short range system. On the other hand it should be possible to observe a "reminiscence" of 
the mean field Td > T c in the 1RSB case by looking at the energy density relaxation when the 
system, starting from a random configuration, is quenched abruptly (i.e. infinite cooling rate) 
to a low T value. 

Recent numerical results on a glass-forming system 



22 seem to show that Td is identi- 



fiable in short-range models as the temperature below which the energy density e(t) displays 
a two step relaxation. The first step corresponds to the relaxation of the system to some 
metastable states with a mean field like behaviour |L2], e(t) oc t~ a where a is weakly depend- 
ing on T. The second one, that happens on a remarkably larger time scale, is due to the slow 
decay of metastable states by activated processes. 

We study the relaxation behaviour of e(t) by considering a non small size L=25. Starting 
from a random initial configuration the system is quenched to the final temperature and it 
evolves by Monte Carlo, each spin being suggested to flip in one MC step. Presented data are 
averaged over two independent dynamics. 




0.2 0.4 0.6 0.E 

.-0.1 7 



0.2 0.4 0.6 0.8 1 0.2 0.4 0.6 0.8 

.-0.15 + -0.1 



Figure 8: Data on e as a function of y = t de /^( T ) at different temperatures. From left to right 
and top to bottom, T=2.6(+) , 2.3(x), 1.95(*) and 1.3(D), where y=t~ - 20 , t~ - 17 , t" - 15 and 
£-0.10 reS p ec tively. 



The energy density e is presented in [Fig. 8] as a function of y—t~" ae ^ T ' at the different 
temperatures considered T=2.60 (~ T c ), 2.30, 1.95 and 1.30. These data are well consistent 
with a linear behaviour of e as a function of y also at large times (small y), i.e. there is no 
evidence for a two step relaxation. Moreover at the higher temperatures we have considered 
the system seems near to the equilibrium at the end of the time window and correspondingly 
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the values of e result compatible with the equilibrium ones obtained from the previously shown 
data on smaller systems. 

The observed behaviour gives evidence for the dynamical critical temperature coinciding 
with the statical one, = T c , as it happens in models with FRSB. 

Since we are looking at a second order transition, a is related to the dynamical critical 
exponent z by a e (T) = d e /z(T), where d e — d — 1/v. We get a e ~ 0.2 at T = 2.6 ~ T c that 



agrees well with the previous estimates [0 v ~ 2/3, z(T c ) ~ 7. Our statistics is inadequate to 
give meaningful results on the behaviour of z(T), which appears however definitely dependent 
on T, the obtained estimation being consistent with z(T) = z(T c ) T c /T. 
To conclude this section we present in [Fig. 9] data on 

i N i M 

^(2M) = ^E<T7E^(2ty/(t)> (22) 
iV i=i m i=i 

plotted as a function of y = t~~ aq c, where a Qc = a q (T ~ T c ) ~ 0.06. 

At T = 2.6 ~ T c the observed behaviour is consistent with C going to zero linearly in 
t~ aq c. By scaling laws one has a qc = d q /z(T c ), where d q = (d — 2 + r])/2. We note that 
a Qc ~ 0.06 agrees with 77 ~ 0, z(T c ) ~ 7. Our statistics results inadequate to fit data at 
lower temperatures where the behaviour results consistent with C going to a non-zero value 
for t — >• 00 at T < T c . 
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Figure 9: Data on C(2t, t) as a function of y = r d «/* (Tc) = t" - 06 at T=2.6(+) , 2.3(x), 1.95(^ 
and 1.3(D). 
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5.2 Aging and generalized FDT 



The data presented in this section are obtained by MC simulations on a L=16 system at 
T—2.0, i.e. well below the critical temperature. We average over 50 different dynamics (initial 
conditions and realizations of the thermal noise). Starting from a random configuration, 
the system is quenched abruptly to the final temperature. After t w steps we apply a small 
perturbing magnetic field 8h = h = 0.4, measuring the two times correlation function 

^ m 

C(t + t w , t w ) = (—J2 + t w )<Ti l (t w )) (23) 
M i=\ 

and the staggered magnetization 

■y N -y M 

m s [h]{t + t w ) = - £<— £ a\(t + t w )). (24) 
iv i=1 ivi l=1 

The behaviour of C(t + t w ,t w ) for different values of t w is shown in [Fig. 10] on the left. 
On the right we plot the same data as function of t/t w . The hypothesis of simple aging, i.e. 
C(t + t w ,t w ) = C{t/t w ), seems to be not well verified. 
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Figure 10: Data on C(t + t w , t w ) plotted as function of t (left) and as function of t/t w (right) 
for t w =2%+), 2 8 (x), 2 10 (*), 2 n (a), 2 12 (A), 2 13 (Q), 2 14 (.). 



In [Fig. 11] we plot m s [h](t + t w )/(/3 ho) as a function of C(t + t w ,t w ) for different values 
of t w . Here the expected scaling is quite well satisfied and we are therefore allowed to take 
the large times behaviour of this quantity as an estimation for S(C). The behaviour of 
m s [h](t + t w ,t w )/(P ho) at t w = 2 17 (i.e. the largest value considered) is presented in [Fig. 
12] together with the estimation of y(q) = J 1 dq' Jq dq"P(q") as obtained by the L=5 data on 
P(q). As expected, the dynamical function S(C) results perfectly compatible with the statical 
one y(q). This dynamical quantity, obtained by simulating a non small system, gives further 
evidence that the model does not behave like MF 1RSB spin glasses. 

We have shown in the figure the line corresponding to 1 — C just to emphasize the small 
t region where FDT is satisfied. The qi value is evaluable as the one where y begins to move 
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Figure 11: Data on m s [h](t + t w )/((3 h ) plotted as function of C(t + t w ,t w ) for t w =2 e (+), 
2 8 (x), 2 10 (*), 2 11 (D), 2 12 (A), 2 13 (Q), 2 14 (.). 




C,q 



Figure 12: Data obtained from the dynamics on m s [h](t + t w )/(f3 ho) plotted as function of 
C(t + t w ,t w ) at t w =2 17 (+) compared with the behaviour of y(q) as obtained from the P(q) for 
L=5 (the line with error bars). The * represents a dynamical estimation for lim^o S(C) as 
obtained by a simulation in which the magnetic field is switched on from the beginning. We 
also plot 1-C (see the text). 



17 



away from this line. In mean field 1RSB models we observe a straight line for y going from qi 
to (the dynamical function S(C) showing the same structure than the statical one): 

„ ( „\ _ / (! -™)(<?i -q) + (i -qi) for q< qi , . 

yW-\ l- q for q> Ql [Zb) 

This behaviour has been recently observed in simulations of structural glasses 0, by defining 
appropriately the correlation function C and the response function to small perturbations R. 

In the case of our model, instead of been in agreement with (p5j), S(C) (we dynamically 
estimate 5(0) extrapolating data obtained by a simulation in which h has been switched on 
from the beginning) results more similar to the one observable in FRSB models [|T7| [23|, 
confirming the results on the equilibrium probability distribution P(q) at low temperatures. 



6 Discussion and conclusions 

Summing up our main results: 

• The transition appears to be of second order, characterized by a diverging correlation 
length for T — > T c ~ 2.6. Our estimations for the critical exponents agree with the ones 



previously obtained [12[ v ~ 2/3, r\ ~ 0, z ~ 7. 



• There is no evidence for the dynamical critical temperature Td being greater then the 
statical one, where in finite dimensional models that are 1RSB in mean field Td > T c 
should mark the onset of two steps relaxation processes, . 

• The distribution probability of the overlaps P(q) results non trivial at low temperature 
but its behaviour is very different from the one of mean field 1RSB models. 

• Well below T c , we find a P(q) still definitely non-zero on the whole interval [— q\,q-\\- 
This result is confirmed by out of equilibrium simulations on a quite larger lattice (i.e. 
by the behaviour of S(C) which is related to P(q) by the generalized FDT relation). 

• The behaviour of P(q) seems nevertheless also not compatible with the one usually 
encountered in FRSB models, this resulting well evident if one looks at the cumulant 

g(T). 

Our possible interpretation of this phenomenology is the following. In this model the 
disorder is not space- invariant. In each particular sample there will be a finite density of 
regions which are more likely than the rest of the system to freeze into a SG phase. Hence 
when T > T c there is a non- zero probability to find a region of space where the system is in 
the 1RSB phase. The typical size of these regions will diverge when T approaches T c since 
the system freezes at the critical temperature. In these regions the spins are very correlated 
and one can formally divide the total high temperature spin glass susceptibility into the 
contribution of bubbles of frozen spins and the contribution of the rest of the system. We 
write 

Xsg = Xb + Xr (26) 
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where B stands for "bubbles" and R for "rest" . So, for T approaching T c from above one has 

Xb^(T-T c )-^ XR = 0(1) (27) 

It is important to note that this kind of transition, at least from this analytical point of 
view, is very different from an ordinary second order transition in mean field, where a zero 
mass mode on the mean field solution causes diverging correlations. Here the high temperature 
solution is stable and non-perturbative effects dominate the approach to the transition. Within 
the bubbles the transition can still be similar to the MF 1RSB. It is the weight of the bubbles 
that grows continuously and diverges when T — ► T c . 

It is possible that a similar phenomenon is responsible for the fact that the function P(q) 
does not tend to the sum of two delta functions at T\T C (opposite to what happens in mean 
field theory). This point deserves a more careful study. 
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